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It has been recently proved that a slab of negative refractive index material acts as a perfect lens in 
that it makes accessible the sub-wavelength image information contained in the evanescent modes 
of a source. Here we elaborate on perfect lens solutions to spherical shells of negative refractive 
material where magnification of the near-field images becomes possible. The negative refractive 
materials then need to be spatially dispersive with e{r) ^ 1/r and ^(r) ~ 1/r. We concentrate on 
lens-like solutions for the extreme near-field limit. Then the conditions for the TM and TE polarized 
modes become independent of ^ and e respectively. 

I. INTRODUCTION 

The possibility of a perfect lens whose resolution is not limited by the classical diffraction limit has been subject 
to intense debate by the scientific community during the past two years. This perfect lens could be realised by using 
a slab of material with e = /i = — 1 where e is the dielectric constant and the fx is the magnetic permeability. Veselago 
had observed!^ that such a material would have a negative refractive index of n = —^fefi = — 1 (the negative sign of 
the square root needs to be chosen by requirements of causality), and a slab of such a material would act as a lens in 
that it would refocus incident rays from a point source on one side into a point on the other side of the slab (See Fig. 
1). Due to the unavailability of materials with simultaneously negative e and /i, negative refractive index remained an 
academic curiosity until recently when it became possible to fabricate structured meta-materials that have negative e 
and /X 0,0,0. Most of the negative refractive materials(NRM), so far, consist of interleaving arrays of thin metallic 
wires (that provide negative e 6]) and metallic split-ring resonators (that provide negative \i 0). Although some 
initial concerns were expresse d that the observed effects in these experiments were dominated by absorption, the 
recent experiments of Refs. [9I ITol illl W \ have confirmed that negative refractive materials are today's reality. 

It was demonstrated by one of us that the NRM slab acts a lens not only for the propagating waves (for which 
the ray analysis of Veselago is valid) but also for the evanescent near-field radiation . This phenomenon of perfect 
lensing becomes possible due to the surface plasmon states lial that reside on the surfaces of the NRM slab which 
restore the amplitudes of the decaying evanescent waves 
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FIG. 1: Radiation from a point source on one side of a slab of material with i e = —1 and fj, = —1 is refocussed into a point on 
the other side. The rays representing propagating waves are bent on to the other side of the normal at the interfaces due to 
the negative refractive index of the slab. 
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numerical (FDTD) simulations that an incident pulse is temporarily trapped at the interfaces for a considerable time 
[T9l |. For a detailed description of the perfect slab lens, we refer the reader to Ref. P,0,|23|- The 'perfectness' of 
the perfect lens is limited only by the extent to which the constituent NRM are perfect with the specified material 
parameters. Absorption in the NRM and deviations of the material parameters from the resonant surface plasmon 
conditions of the perfect lens causes significant degradation of the subwavelength resolution possible 14, 21, 22, 23j]. 
We have sug gested some possible measures to ameliorate this degradation of the lens resolution by stratifying the lens 
medium |24| and introducing optical gain into the system f2j |. 

The image formed by the NRM slab lens is identical to the object and hence there is no magnification in the 
image. Lenses are mostly used to produce magnified or demagnified images and the lack of any magnification is a 
great restriction on the slab lens on which most of the attention in the literature has been focussed. The slab lens is 
invariant in the transverse directions and conserves the parallel component of the wave- vector. To cause magnification 
this tranverse invariance will have to be broken and curved surfaces necessarily have to be involved. The perfect lens 
effect is dependent on the near-degeneracy of the surface plasmon resonances to amplify the near-field, and curved 
surfaces in general have completely different surface plasmon spectrum It was recently pointed out by us that 
a family of near-field lenses (in the quasi-static approximation) in two-dimensions can be generated by a conformal 
mapping of the slab lens 20]. Thus a cylinderical annulus with dielectric constant e = —1 was shown to have a 
lens-like property of projecting in and out images of charge distributions. Similiarly in Ref. and |2^, it was 
shown how a general method of co-ordinate transformations could be used to map the perfect slab lens solution for 
the Maxwell's equations into a variety of situations including the cylinderical and spherical geometries respectively. 

In this paper, we elaborate on the perfect lens solutions in the spherical geometry and show that media with 
spatially dispersive dielectric constant e{r)/sim\/r and magnetic permeability iJ,{r) / siml / r can be used to fabricate 
a spherical perfect lens that can magnify the near-field images as well. In section-2 of this paper, we will present 
these perfect lens solutions of the Maxwell's equations for the spherical geometry. In section-3, we will examine the 
solutions in the extreme near-field limit or the quasi-static limit which is useful when the lengthscales in the problem 
are all much smaller than a wavelength. Then the requirements for TM and TE polarizations depend only on e l/r^ 
or ~ 1 /r^ respectively. This is useful at frequencies where we are able to generate structures with only one of e or 
fi negative. We will investigate the effects of dissipation in the NRM and point out the connections to the 1-D slab 
lens solutions. We will present our concluding remarks in Section-4. 



II. A PERFECT SPHERICAL LENS 



Consider a spherically symmetric system shown in Fig. 2 consisting of a spherical shell of NRM with the dielectric 
constant S-(r) and fJ—(r) imbedded in a positive refractive material with e+(r) and fj,-^-{r). First of all we will find 
the general solutions to the field equations with spatially inhomogeneous material parameters: 

V X E = «w/ioAi(r)H, V x H icjeo£(r)E (1) 

V • D = 0, V • H 0, (2) 
D = e(r)E, B = ^(r)H. (3) 

Under these circumstances of spherical symmetry, it is sufficient to specify the quantities (r • E) and (r • H) which 
will constitute a full solution to the problem. Let us now look at the TM polarised modes r • H = 0, implying that 
only the electric fields have a radial component Er. Operating on Eqn. (1) by V, we have 



V X V X E iLu^oV X [Ai(r)H] , 

/ X / ^^ V//(r) 
= ^Mr)s(r)E + ..-i^^ 



i^/Li(r)e(r)E + iio-^ x V x E. (4) 



Using Eqns.(2) and (3) we have 

V • D = V • [e(r)E] = Ve(r) • E + £(r) V • E = 0, (5) 

and if we assume e(r) — e{r) and /i(r) = tJ-{r), we have 

e'(r) e'(r) 
re(r) re[r) 

We note the following identities for later use: 

V X V X E = V(V • E) - V^E, (7) 
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FIG. 2: A spherical shell with negative e-{r) — 1/r and fJ.-(r) ~ — 1/r images a source located inside the shell into the 
external region. The media outside have postive refractive index, but e_(r)i ~ 1/r and ^J.-{r) ~ 1/r. The amplification inside 
the spherical shell of the otherwise decaying field is schematically shown. 



and using Eqn. (6) we also note that 



r-V(V-E) r-V 



V2(r-E) = r- V2E + 2V-E, 

e'ir) 



(8) 



e{r) 
d fe'ir) 



-Er 



Or \ e{r) 
d fe'ir) 



(rEr 



e'ir) 



Er 



(9) 



dr \ e(r) / \ e(r) 

We now take a dot product of r with Eqn. (4), and use the Eqns. (6), (7), (8) and (9) to get an equation for (riJ^) as: 



V^irEr)- 



d_ 

dr 



e'ir) 
e{r) 



irE,.) 



+ ^irEr)+eir)f,ir)^irEr)=0. 
reir) 



(10) 



This equation is separable and the spherical harmonics are a solution to the angular part. Hence the solution is 
irEr) = U{r)Yimid, (j>) where the radial part t/(r) satisfies 



1 d f^2dU^^lil + l)^ ^ d 



e'ir) 



eir) 



U 



e'ir) iJ^ 
+ ^C/ + e(r)Mr)-C/ = 0. 
reir) cJ- 



^2 Qr \ dr ) ' dr 

If we choose e(r) — ar^ and /i(r) = /5r'', we can have a solution U (r) r" and we get 

[n(n +l)-lil + l)+pin-l)+ p]r"-^ + a(3uj^ /c\p+''+'' = 0, 

implying p + q — —2 and 



n± = l/2 -ip + 1) ± y/ip + 1)2 + nil + 1) - ia(3ujyc 
Hence the general solution can be written as 

Erir) = J2 [n+Ai^r^-^-^+n.Bi^r^—^] y,™(0,0) 

l.m 



(11) 



(12) 



(13) 



(14) 
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and 



H(r) = J2 + Bimr''-]Xim{e, Cj)), 



where the vector spherical harmonic 



(15) 



(16) 



A simihar solution can be obtained for the TE modes with r • E = 0. 

Now assuming an arbitrary source at r = oo, we can now write down the electric fields of the TM modes in the 
different regions for the negative spherical shell of Fig. 2 as 



EW(r) = 5^[n+4;2^"+-i+n_B, 



(1), 

Im' 



l.m 



E(3)(r) = E[n+4^V"+-i+n_B, 



Im ' 



ao < r < ai, 
ai < r < 02, 
02 < r < oo, 



(17) 
(18) 
(19) 



l.m 



and similiarly for the magnetic fields. Note that the correspond to the field components of the source located 

— 0. Now the tangential components of the magnetic fields and the normal 
have to be continuous across the interfaces. Under the conditions » = — 1, 



at r = flo. For causal solutions A 
components of the displacement fields 
q = -1, e+(ai) = -e-(ai) and £+(02) 



= — e_ 


(02), 




= 0, 


.(2) 


-{ 


Im 


= 0, 


bP 

Im 


-{ 



V/(i+l)-a/3a;2/c2 



B 



(1) 



2^^l{l+l)-oil3uj^/c 



Im 

The lens-like property of the system becomes clear by writing the field outside the spherical shell as 

£;(3) = : 



B\2Yim{e,<i>). 



(20) 

(21) 

(22) 

(23) 
(24) 



Hence apart from a scaling factor of 1/r, the fields on the sphere r = 03 = (a|/af)ao are identical to the fields on the 
sphere r = ao- There is also a spatial magnification in the image by a factor of 02/ a\. 

Let us note a couple of points about the above perfect lens solutions in the spherical geometry. First, for r > 0,3. i.e. 
points outside the image surface the fields appear as if the source were located on the spherical image surface (r — 03). 
However, this is not true for points a2 < r < within the image surface. Second, given that £-(02) = —£+(02), 
we have the perfect lens solutions if and only if n+ = — n_ which implies that p = —1 in Eqn. (13). Although 
the solutions given by Eqn. (14) occur in any medium with en ~ l/r^, the perfect lens solutions only occurs for 
£ ~ ^ ~ 1/r. Here we have written down the solutions for the TM modes. The solutions for the TE modes can be 
similiarly obtained. 



III. THE SPHERICAL NEAR-FIELD LENS 



As it has been pointed out in the previous section, the 'power' solutions arc good for any s{r) ^ and /i(r) ^ 
such that p + q = —2. However the perfect lens solutions for the Maxwell's equations result only for the single case 
of p = q = —1. In the quasi-static limit of w ^ and / » we can relax this condition. In particular, by 

setting e{r) ~ 1/r^ and ^{r) = constant, we can have a perfect lens for the TM modes alone. Similiarly, we can have 
a perfect lens for the TE polarization by having /z ~ 1/r^ and e = constant. 
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This extreme near-field limit is both important and valid for situations when all lengthscales in the problem are 
much smaller than a wavelength of the radiation. This becomes useful at frequencies where we can only generate 
media with either negative e and positive /i, or, negative fj, and positive e. Examples are the silver slab lens at optical 
frequencies [ij, the imeta-materials (Swiss rolls) used for MRI at radio- frequencies '2^. Particularly at radio- and 
microwave frequencies, we currently can practically engineer the required meta-matcrials with spatially dispersive 
characteristics at the corresponding length scales. Further it also lifts the restriction that the system has to have a 
spatially dispersive material parameters even outside the spherical shell of NRM. In this section we will work in this 
extreme near field limit. Then it is sufficient to solve the Laplace equation and we present lens- like solutions to the 
Laplace equation below. 

Consider the spherical shell in Fig. 2 to be filled with a material with €2(7") ~ — C/r^ with the inner and outer 
regions filled with constant dielectrics of ei and £3 respectively. Let — I everywhere. Now place a charge +q at the 
centre of the concentric spheres and a charge —q at a distance ag from the centre inside region-1. We will consider 
the z-axis to be along the dipole axis and make use of the azimuthal symmetry here, although it is clear that our 
results do not depend on any such assumption of azimuthal symmetry. Thus all our charge and their images will now 
lie along the Z-axis. 

Now we will calculate the potentials in the three regions that satisfies the Laplace equation and the continuity 
conditions at the interfaces. The potential in region-l(r < ai) can be calculated to be (using the azimuthal symmetry): 



Vi(r) 



47rei 



E 

1=1 



Aur^Pi{cosl 



Y Pi {cost 



(25) 



Note that the second term in the above expansion arises due to the dipole within the sphere. It can be shown (see 
appendix-1), that the general form of the potential in region-2 (ai < r < 02), where the dielectric constant varies as 
l/r^ is 



" 1=1 

In region-3(r > 02), the potential is given by: 

^3(r) = 



47re3 



E 

1=1 



(26) 



(27) 



Now we must match the potentials at the interfaces at r = ai and r = 02 (Put eo = 1) to determine the A and B 
coefficients. The conditions of continuity of the potential and the normal component of D at the interfaces are 



l/i(ai) = 1/2(01) 

dVijai) dV2{ai) 
ei 7; = £2 7; 



dr dr 
We determine the coefficients from these conditions to be: 



^2(02) = V'3(a2), 
9^2(02) dVsia^) 



£2- 



dr 



£3- 



dr 



(28) 
(29) 



All = 
A21 = 
B21 = 



B: 



31 



a- 


h I)a[,{[/e2(a2) ~ [1 + l)e3][ei + £2(01)] - 


[e2(a2) + e3][(? + l)ei 


-/e2( 


„2i + l 

ai)] Ji+i} 
"1 


;(H 


-l)[£i 


+ e2(ai)][e2(a2) + esla^'^^ + [^ei - 
{2l + l)[le2{a2)-{l 


[1 + I)€2{ai)][le2{a2) 
+ l)e3]aoar^ 




iN 1 2i-|-l 


1(1^ 




+ e2(ai)][e2(a2) + £3102'+^ + [^ei - 
(2/ + l)(/-|-l)[e2(a2)^ 


[1 + I)€2{ai)][le2{a2) 
- e3ja2aoai 


-a + 


iN 1 2i + l 


1(1^ 




+ e2(ai)][e2(a2) + £3102'+^ + [^ei - 
(2/ + l)2e3e2(a2)a 


{1 + I)e2{ai)][le2{a2) 

„2(i + l) -1 
0"2 "1 


-a + 


l)e3]af+^ 






+ e2(ai)][e2(a2) + £3102'^^ + [l^i - 


{I + I)€2{ai)][le2{a2) 


-a + 


l)e3]af+^ 



(30) 
(31) 
(32) 
(33) 



Under the perfect lens conditions 



£2(01) = -ei, and £2(02) = -£3, 



(34) 
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we have: 



Au 


= 0, 




A21 


1 


ah 




„2(i+l) ' 

"1 


B21 


= 0, 




B31 


_ ^3 
Si 


/„ \2(;+i) 



Hence the potential outside the spherical shell for r > 02 is 



Aires 



£1 



2(1+1) 



which is the potential of a dipole with the positive charge at the origin and the negative charge at 03, where 

2 



03 



and of strength 



52 



£3 

£1 



ao, 



(35) 
(36) 

(37) 
(38) 

(39) 

(40) 
(41) 



as ^a/ei = (01/02)^. Thus, on one side of the image (the region r > 03) the fields of a point charge located at 03 are 
reproduced. However it should be pointed out that there is no physical charge in the image location and, the fields 
on the other side of the image (i.e. in the region a2 <r < a^) do not converge to the fields of the object and cannot 
do so in the absence of a charge in the image. Further there is no change in the strength of the charge either. There 
is a magnification in the image formed by a factor of (02/01)^. 

Now let us consider the case of a point source placed at 03 in the outer region. Again assuming the z-axis to pass 
through as, we can write the potentials in the three regions as 



Viir) 
V2{v) 
Vs{r) 



+q 

A'K£i 



Y,Aur'Pi{cose) 



V 



r < ai, 



+q 

47r£3 



1=0 

00 



E 

1=0 

00 



A2ir 



i+i 



B21 



E 



B 



31 



J+1 



P((C0S( 

Pi (cos 6 



V 



ai < r < 02, 



a2 <r < as, 



(42) 
(43) 

(44) 



where the first term in Vs{r) comes from the point source at 03. Now applying the conditions of continuity at the 
interfaces, we can similiarly obtain for the coefficients as before. In the limiting case of £2(^1) = —£i and £2(^2) = —£3, 
we have 



Hence the potential inside the inner sphere is 



Vi{v) = 



All 
A21 
B21 
Bsi 

q 

A'K£\ 



£1 
£3 

0, 

1 



1 



£3 a- 
0. 



21+1 ■ 



E 



£1 
£3 



Pi (cos 61), 



(45) 
(46) 
(47) 
(48) 

(49) 



i.e., that of a point charge of strength q\ = {£1/ £s){0'i/ 0-2)'^ = g at ao = 03(01/02)^. As before, for the inner region 

of r < oq, the system behaves as if there were a single charge of strength q located at r = ao- Thus the shell has a 
lens-like action. We note that there is a demagnification of (ai/a2)^ in this case. 
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A. Similiarities to the 1-D slab lens 



Let us point out the similiarities to the planar slab lens. In both cases, the electromagnetic field grows in 
amplitude across the negative medium when the perfect lens conditions are satisfied at the interfaces: as an 
exponential(exp[+/ca;2:]) in the planar lens and as a power of the radial distance r' in the spherical lens. The de- 
caying solution away from the source is absent in the negative medium in both cases. Further, when the perfect 
lens conditions are matched at both the interfaces, there is no reflected wave in both the planar slab as well as the 
spherical lens: i.e. the impedance matching is perfect as well. In addition this mapping preserves the strength of the 
charge. 

The key differences, however, are the different dielectric constants on either sides of the spherical shell of the negative 
medium. This is a direct consequence of the spatial dependence of the negative dielectric constant which relates 
the two positive dielectric constants to be ei — (ai/a2)^e3. But this need not be a particular restriction as we can use 
the ideas of the asymmetric lens to terminate the different positive media at some radii bevondfl^. The net result is 
that the image can now be magnified (or demagnified) when the image of the charge (source) is projected out of (or 
into) the spherical shell, which is true in the 2D cylinderical lens as well poj |. 



B. Possibility of the asymmetric lens 

In the case of a planar slab, it was possible to have the perfect lens effect by satisfying the required conditions at 
any one interface - not necessarily at both interfaces p7|. Particularly , in the limit of very large parallel wave vectors 
the lensing is indeed perfect, although the image intensity differed from the source by a constant factor. Similiarly 
let us now investigate the effects of having the perfect lens conditions in the case of the spherical lens at only one of 
the interfaces. 

Let us consider first, the case of projecting out the image of a point source from inside the spherical shell to outside 
and enable the perfect lens conditions only at the outer interface €2(7" = 02) = —£3 and have an arbitrary ci. Now 
the A and B coefficients come out to be 

af+Mei- (/ + l)£2(ai) ^ ' 

_ (2/ + 1)4 

- 7- „ , . ^ 2/4-2 (^^) 





{1 + 


l)[ei +£2(01)] 




lei - 


- (/ + l)£2(ai) 




{21 


+ 1)4 


lei — 


(1 + 


l)e2{ai)af+' 











(2Z4 


- 1)£34 f 



B21 = (52) 

'a \ 

" lei{ai/a2r+ll + l)es ['^J ^^^^ 

Only the growing solution within the negative spherical shell remains. The coefficient of the decaying solution {B21) 
remains strictly zero. Thus, amplification of the decaying field at least is possible in this case as well. But there is a 
finite reflectivity in this case. However, the solution outside for r > a2 is not the exact image field of the point source 
as the coefficient B31 has an extra dependence on / through the dependence on the dielectric constants. Moreover, the 
process does not preserve the strength of the charge due to the different dielectric constants involved. This should 
be compared to the solution of the planar asymmetric slab lens where, at least in the electrostastic limit, the system 
behaved as a perfect lens. In this case, the system behaves as a spherical asymmetric perfect lens only in the limit 
of large Z — > cx). The solution outside the spherical shell is the same when we meet the perfect lens condition on the 
inner interface - just as in the case of the planar slab lens. However, the reflection coefficient is again non-zero, but 
different to the earlier case. In either case, the fields are largest at the interface where one meets the perfect lens 
conditions or the interface on which the surface plasmons are excited. 



C. Effects of dissipation 



Media with negative real part of the dielectric constant are absorptive (as all metals are), and hence we can write 
the dielectric constant e{r) — C/r^ + iei{r) (Note that ^ as well for us to be able to write the solution in 
the following form). Consider the first case of projecting out the image of a dipole located within the spherical shell 
where the potential outside the shell is given by Equation (27) and -83; is given by Equation (33). When we have a 
dissipative negative medium and have the perfect lens conditions at the interfaces on the real parts of the dielectric 
constant alone, £2(01) = — £1 -I- ieiiai) and £2(02) = —£3 + *£i(a2). In parallel with the case of the planar lens, we 
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note that the denominator of -63; consists of two terms, one containing a power of the (smaller) radius ai and the 
other containing a power of the (larger) radius 02- Crucially the amplification of the evanescent fields depends on the 
possibility that the smaller power dominates by making the coefficient of the larger term as close as possible to zero. 
The presence of the imaginary part of the dielectric constant would not allow the coefficient to be zero and the image 
restoration is good only as long as the term containing ai dominates in the denominator of i.e., 

l{l + l)e,iai)e,{a2)al'+^ « [{21 + l)ei - i{l + l)e,(ai)] [-(2^ + 1)^3 + ^e^{a2)]af+^ (54) 

Hence a useful estimate of the extent of image resolution can be obtained by noting the multipolc I for which the two 
terms in the denominator are approximately equal,14j. We obtain for this value: 

, ^ ln{3£i£3/[£i(ai)£i(a2)]} , . 

'max — o 1 / / \ • vOOj 

2 In (02/01) 

Higher order multipoles are essentially unresolved in the image. We can similiarly obtain the same criterion by 
considering the second case of tranferring the image of a charge located outside the spherical shell into the inner 
region. Again, we can consider the effects of deviating from the perfect lens conditions on the real part of the 
dielectric constant as well and obtain a similiar limit for the image resolution due to those deviations. 



IV. CONCLUSIONS 



In conclusion, we have presented a spherical perfect lens which enables magnification of the near-field images. The 
perfect lens solution requires media with e{r) ~ l/r and fj,{r) ^ 1/r and the conditions £-(ai 2) = —£+(01^2) Sind 
/^-(ai,2) = ~A'+(ai.2) at the interfaces of the spherical shell of the NRM. We have shown that in the quasi-static limit 
of small frequencies {lo — > 0) and high-order multipoles Z 3> \p\, this condition can be relaxed and the two polarizations 
(TE and TM modes) decouple. Thus a shell with negative dielectric constant £-(r) — with fi = constant can 
act as a near-field lens for the TM polarization while fJ.-{r) ^ — l/r^ with £-(r) = constant acts as a near-field lens 
for the TE modes. We have shown that dissipation in the lens material, however, prevents good resolution of higher 
order multipoles. Thus while the near-field lenses work best for the higher order multipoles, dissipation cuts-off the 
higher order multipoles. Further the spherical lens works in the asymmetric mode only in the limit of high order 
multipoles. Thus one has to find an intermediate regime where dissipation does not wipe out the near-field image 
information and yet the meta-materials work. This is the design challenge involving these near-field lenses. 



Appendix: Solution of the Laplace equation in a spatially varying medium 



We have to solve the Maxwell's equations in material media 

V • D = 0, => V(£E) = 
Using E = — VU, where V{r) is the electrostatic potential we have: 

£(r)V2U-|-V£(r)-VU = 



(56) 



(57) 



If £(r) has only a radial dependence (as in our case ~ ^/r'^), then V£(r) — f{de/dr) and we can separate the solution 



as U(r) = 



where the Y/m is the spherical harmonic and the radial part U (r) satisfies: 



£(r) 



<PU_ l{l + l) 



U 



de 
dr 



dU 
dr 



= 



(58) 



To have a solution as a single power of r, the only choices possible for the dielectric constant are either £ = C, a 



constant - the usual case , or £ = C jr . In the latter case the solution is U{r) 
solution can then be written as 



or \J(r) - r-C-i). The full 



U(r) = [^'"^'^' 



(59) 



i=0 



[1] J.B. Pendry, Phys. Rev. Lett. 85, 3966 (2000). 



9 



[2] V.G. Veselago, Usp. Fiz. Nauk 92, 517 (1964) [Sov. Phys. Uspekhi, 10, 509 (1968).] 

[3] D.R. Smith, W.J. Padilla, D.C. Vier, S.C. Nemat-Nasser and S. Schultz, Phys. Rev. Lett., 84, 4184 (2000). 
[4] R.A. Shelby, D.R. Smith and S. Schultz, Science, 292, 77 (2001). 
[5] A. Grbic and G.V. Eleftheriades, J. Appl. Phys. 92, 5930 (2002). 

[6] J.B. Pendry, A.J. Holden, W.J. Stewart and I. Youngs, Phys. Rev. Lett., 76, 4773 (1996); J.B. Pendry, A.J. Holden, D.J. 

Robbins and W.J. Steward, J. Phys.: Condens. Matter, 10, 4785 (1998). 
[7] J.B. Pendry, A.J. Holden, D.J. Robbins and W.J. Steward, IEEE trans. MTT, 47, 2075 (1999). 
[8] N. Garcia and M. Nieto-Vesperinas, Opt. Lett. 27, 885 (2002). 

[9] C.G. Parazzoh, R.B. Greegor, K. Li, B.E.C. Kotenbah and M. Tanielan, Phys. Rev. Lett. 90, 107401 (2003). 
[10] A.A. Houck, J.B. Brock, LL. Chuang, Phys. Rev. Lett. 90, 137401 (2003). 
[11] A. Grbic and G.V. Eleftheriades, AppL Phys. Lett. 82, 1815 (2003). 
[12] P.V. Parimi, W.T. Lu, P. Vodo, J. Sokoloff and S. Sridhar, cond-mat/0306T09l 
[13] H. Raether, Surface Plasmons, (Springer- Verlag, Berlin, 1988) 

[14] D.R. Smith, D. Schurig, M. Rosenbluth, S. Schuhz, S.A. Ramakrishna and J.B. Pendry, Appl. Phys. Lett. 82, 1506 (2003). 

[15] F.D.M. Haldane, cond-mat/0206420 

[16] G. Gomez-Santos, Phys. Rev. Lett., 90, 077401 (2003). 

[17] S.A. Ramakrishna, J.B. Pendry, D.R. Smith, D. Schurig and S. Schuhz, J. Mod. Optics, 49,1747 (2002) 
[18] X.S. Rao and C.N. Ong, Phys. Rev. B 68, 113103 (2003). 

[19] S. Foteinopolou, E.N. Economou and CM. Soukoulis, Phys. Rev. Lett. 90 107402 (2003). 

[20] J.B. Pendry and S.A. Ramakrishna, J. Phys.: Condens. Matter 14, 8463 (2002). 

[21] J.T. Shen and P. Platzmann, Appl. Phys. Lett. 90, 3286 (2002) 

[22] Z. Ye, Phys. Rev. B 67, 193106 (2003). 

[23] N. Fang and X. Zhang, Appl. Phys. Lett. 82, 161 (2003). 

[24] S.A. Ramakrishna, J.B. Pendry, M.C.K. Wihshire and W.J. Stewart, J. Mod. Optics 50 1419 (2003). 
[25] S. A. Ramakrishna and J.B. Pendry, Phys. Rev. B. 67, 201101(R) (2003). 
[26] V.V. Klimov, Opt. Commun. 211, 183 (2002). 
[27] J.B. Pendry, Opt. Express 11, 755 (2003). 

[28] J.B. Pendry and S.A. Ramakrishna, J. Phys.: Condens. Matter 15, 6345 (2003) 

[29] M.C.K. Wiltshire, J.B. Pendry, l.R. Young, D.J. Larkman, D.J. Gilderdale and J.V. Hajnal, Science 291, 848 (2001); Opt. 
Express, 11, 709 (2003). 



